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ABSTRACT 



> 

, X-ray line profile analysis has proved to be the most direct diagnostic of the kine- 

1/^ ■ matics and spatial distribution of the very hot plasma around O stars. The Doppler- 

O^l ' broadened line profiles provide information about the velocity distribution of the hot 

, plasma, while the wavelength-dependent attenuation across a line profile provides in- 

formation about the absorption to the hot plasma, thus providing a strong constraint 
on its physical location. In this paper we apply several analysis techniques to the 
r"| . emission lines in the Chandra HETGS spectrum of the late-0 supergiant C Ori (09.7 

1' lb), including the fitting of a simple line-profile model. We show that there is distinct 

evidence for blue shifts and profile asymmetry, as well as broadening in the X-ray 
emission lines of C Ori. These are the observational hallmarks of a wind-shock X-ray 
source, and the results for ^ Ori are very similar to those for the earlier O star, ( Pup, 
^ , which we have previously shown to be well-fit by the same wind-shock line-profile 

model. The more subtle effects on the line-profile morphologies in C Ori, as compared 
to C Pup, are consistent with the somewhat lower density wind in this later O super- 
, giant. In both stars, the wind optical depths required to explain the mildly asymmetric 

^ • X-ray line profiles imply reductions in the effective opacity of nearly an order of mag- 

I nitude, which may be explained by some combination of mass-loss rate reduction and 

large-scale clumping, with its associated porosity-based effects on radiation transfer. 
In the context of the recent reanalysis of the helium-like line intensity ratios in both 
C Ori and C Pup, and also in light of recent work questioning the published mass-loss 
rates in OB stars, these new results indicate that the X-ray emission from ( Ori can be 
understood within the framework of the standard wind-shock scenario for hot stars. 

Key words: line: profiles - stars: early-type - stars: mass loss - stars: winds, outflows 
- stars: individual Ori) - X-rays: stars 



1 INTRODUCTION 

X-ray emission from normal, single OB stars has gen- 
erally been explained in terms of shock heating of 
the massive, highly supersonic, radiation-driven winds 
of these very luminous stars dPallavicini et alj Il98l" 
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i n th e context of hot star winds by |Lu£^^^olomoE 
197C ) and later invest igated by iLucv^^^^i^ 



19801): ILuctI il98j): _ lOwocki fc Rvbickil Jl98. 

Owocki. Castor. & Rvbicki (■l988");'Fcldmcior ct al.' (W97sL 
iFcldm cicr. Puis, & Pauldrach (1997b); Runacrcs & Owockil 
i2002f) . The wind-shock X-rays in this scenario arise 
naturally from non-local radiation transfer in the context 
of the standard "CAK" line-dr iven winds of massive stars 
JCastor. Abbott, fc Kleinlll975^ . 

Other wind-shock models of X-ray production have 
also been propo sed, based on co-ro tating interaction regions 
iMullanI Il984) . driven shocks llMacFarlane fc Cassinelli 
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^8£ ) , and inverse Compton scattering JChen fc Whitd 
19911) . Even in the context of the line- force instability mech- 
anism, there are different scenarios based on the self-excited 
instability ( Owocki et al., 1988..') versus the i nstability seeded 
by per turbations at the base of the wind jPeldmeier et alJ 
lliiZH), and one- dimensional simulations versus two- 
dimensional simulations jPessart fc Owockill2003t) . 

Despite this proliferation of models, very few obser- 
vational constraints could be put on any of these wind- 
shock models until recently. This was at least partly be- 
cause of the very limited data quality of X-ray observa- 
tions before the late 1990s. Furthermore, the idea that 
dynamo-driven coronal mechanisms, similar to what is seen 
on the Sun, might be relevant t o hot-star X-ray produc- 
tion, con t inues t o have adherents llCassinelli fc 01sorJll979l: 
Waldron' 198?, 'Smit h et al.) Il993l: IWaldron fc CassineUi 
200 1; Smith et al.., 200'tfl. Indeed, models of surface mag- 
netic field generation and dynamo mechanisms that do not 
involve envelope convection, and therefore might be appli- 
cable to_jnassive_staxSj_Jm^ been proposed (see, 
e.g., iMuUan fc MacDonaldl l|200^ and references therein). 
Speculation on the applicability of such models to mas- 
sive stars, including specifically C, Ori, have, in fact, 
been motivated by the claims of symmetric X-ray emis- 
sion lines and anomalous X-ray line ratios in C Ori 
iMullan fc Waldro n 2006l. Finally, it has recently been sug- 
gested that a hybrid wind-magnetic shock-heating X-ray 
productio n mechanism is in operation on at least som e 
hot stars llGag ne et a l.lll997l: iBabel fc Montmerlelll997allbl: 
l ud-Doula fcO wocki 2 0021 : ISchulz et alj 120031: iGagne et all 

The launch of the Chandra and XMM-Newton tele- 
scopes in 1999, with their high-resolution grating spectrom- 
eters, vastly improved the quality of X-ray spectra available 
from OB stars. Although these missions provided a huge in- 
crease in the amount of information in the X-ray data, they 
have not led to a consensus in the community regarding 
the actual X-ray production mechanism in hot stars. This 
is partly because of the diverse behavior seen in the half- 
dozen or so normal (not interacting binary) O stars thus far 
observed. And it is also partly due to the lack of connection 
between the diagnostics usually employed in the analysis 
of grating spectra of O stars and any quantitative physical 
model. 

The high resolution of the new X-ray grating spectrom- 
eters provides a powerful diagnostic of plasma kinematics 
and location (via the effects of continuum absorption across 
a line), and thus potentially a discriminant among the var- 
ious physical models, in the form of resolved emission line 
profiles. Hot plasma embedded in a fast stellar wind pro- 
duces Doppler-broadened emission lines, with the velocity 
and density structure dictating the detailed form of these 
profiles. Continuum absorption by the unshocked wind is 
stronger on the red sides of emission lines, as the red-shifted 
photons originating in the far side of the wind traverse a 
larger column of material than those emitted from the front, 
blue-shifted side. Overall, then, a wind-shock mechanism, in 
which the shock-heated plasma is embedded in a more-or- 
less spherically symmetric cold wind, should lead to broad- 
ened and asymmetric lines with a blue-shifted centroid and 
a characteristic shape iMacFarlane et alJll99ll : llgnacell200lt 
lOwocki fc Cohenll200ll) . 



The diverse behavior observed in the first hot stars ob- 
served with Chandra and XMM-Newton includes line pro- 
files that are broad, shifted, and skewed in the earliest O 
stars (.Kahn et, al. 2001: Cassinell^^^. 2001') but which are 



quite narrow in early B stars^^^^^^^j200djMew^^l 

3; I; 



j2003l) and young O stars JSchul^^d J l2003in Gagne et al 



2005), with the X-ray emission lines of late O supergiants, 
inclu ding C Ori, having a more intermediate morphology 
( W aldron fc Cassinelljl200ll: iMiller et alJl2002^ . 

The 04 star C, Pup seems to be generally accepted 
as fitting the wind-shock paradigm, based on its line pro- 
files. The broad, shifted, and asymmetric profiles are qual- 
itatively what is expected from a spherically symmetric 
win d source ij CassinBll^^lJ 20Wj). Quantitative analy- 
sis (lKramer~Cohen^^^wocki 12003) shows that the hot 
plasma is distributed throughout the wind above some min- 
imum radius of emission that is approximately half a stellar 
radius above the photosphere; that it is distributed roughly 
as the density-squared of the bulk wind; and that the kine- 
matics of the hot pl asma are consistent with t he underly- 
ing beta- velocity law iLamers fc Cassinellilll99'9l) of the bulk 
wind. This same analysis does, however, show that there is 
significantly less continuum absorption than would be ex- 
pected for a smooth, spherically symmetric wind having 
a mass-loss rate consistent with UV and Hq observations 
and wind opacity consistent with models. This might be 
explained by a reduction in the mass-loss rate or by inho- 
mogeneities in the wind ( "clumping" or "porosity" ) . To the 
extent that the wind-shock picture is applicable to Pup, it 
has generally been supposed, however, that this very early O 
star is the only hot star for which the standard wind-shock 
scenario can explain the Chandra observations. 

The subject of this study, the late-0 supergiant C, Ori, 
meanwhile, has X-ray lines that are broad enough to 
be understood in the conte xt of the wind-shock scenario 
iWaldron fc Cassinellil 1200 J) . It was originally reported, 
however, that there was no systematic trend in the Doppler 
shifts of the emission lines observed with Chandra, and addi- 
tionally, that the forbidden-to-intercombination line ratio of 
Si XIII indicates a location so close to the photosphere that 
it c ould not be explained in the context of wind-shock mod- 
els dWaldron fc CassineUil I2OO1I) . However, no quantitative 
assessment has yet been made of the line profile shapes. In 
this paper we quantitatively examine the shift and asymme- 
try in the X-ray emission lines on C, Ori. We do this first by 
fitting Gaussians to the strong emission lines in the Chandra 
spectrum, and then by performing a non-parametric analy- 
sis of the line shift and asymmetry, and finally by applying 
the simple line-profile model that was successfully used to fit 
the X-ray emission lines in the Chandra spectrum of C, Pup. 
We show that the X-ray emission lines in C, Ori actually can 
be as well fit by standard wind-shock models as those in 
C, Pup, with a similar finding of lower-than-expected wind 
absorption. We also discuss the results of our line-profile 
analysis of C, Ori in light of a reevaluation of the forbidden- 
to-intercombination line ratios that revises the earlier results 
to show no significant conflict with a wind-shock origin for 
the X-rays (Leutcncggcr ct al. 2006;). 

In §2 we briefly describe the observational data and the 
properties of C, Ori. In §3 we assess the blue shifts and skew- 
ness of the line profiles quantitatively but in a non-model- 
dependent way. In §4 we report on fits of an analytic, spher- 
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ically s ymmetric wind emissio n and absorption line-profile 
model lOwocki fc Collenll2001^ to nine lines in the Chandra 
spectrum of C, Ori. In §5 we discuss the results of the model 
fitting and their interpretation, including how these results 
comport with other X-ray diagnostics, especially the helium- 
like forbidden-to-intercombination line flux ratios and UV 
absorption line mass-loss diagnostics. Finally, we summarize 
our conclusions in §6. 



2 THE Chandra DATA AND STELLAR 
PROPERTIES 

The data analyzed in this paper was obtained during the 
Chandra AOl GO phase, using the ACIS-S/HETGS con- 
figuration, and made with nominal pointing at C, Ori. The 
effective exposure time was 73.87 ks, with the data com- 
prising two Obs. IDs: 610 and 1524, taken on 2000 April 8 
and 2000 April 9, respectively. In the combined data, 11,347 
first-order MEG counts were recorded. The dispersed spec- 
trum is quite soft, as can be seen in Figure and there 
were significantly more counts in the MEG than in the HEG 
spectrum, which had only 2508 total first-order counts. We 
therefore used only the MEG spectrum for the line profile 
analysis in this paper. We performed the standard reduction 
and extraction of the dispersed spectra using the basic grat- 
ing threads and CIAO v3.1 and CALDB v2.28. We checked 
the centroids of strong lines separately in the negative and 
positive first-order spectra and did not see any significant 
systematic shift in the wavelengths of the emission lines be- 
tween the negative and positive sides. We wrote the count 
spectra (-1 and +1 orders) to ascii files, and performed the 
analysis with custom-written codes in IDL and Mathemat- 
ica, except for the initial fitting of Gaussian line profiles, 
which we performed within XSPEC. We then repeated the 
fits of wind-profile models to individual lines using a custom 
written model within XSPEC vll.3.1. For the XSPEC fit- 
ting, we used only Obs. ID 610 (exposure time of 59.63 ks), 
as including the second, much shorter, Obs. ID did not im- 
prove the statistics on the fits significantly. For all the model 
fitting reported on in this paper, we used the C statistic to 
assess the goodness-of-fit and parameter confidence limits, 
as the data in the line wings an d nearby co ntinuum have a 
small number of counts per bin llCashlll979D . We discuss the 
fitting procedure in detail in §4. 

The late O supergiant, ( Ori (Alnitak, HD37742, the 
eastern-most of the Ori on belt stars), has a Hipparcos 
distance of 277I49 P^^ iPerrvman et all Il997l) . It has a 
spectral classification of 9.7 and a luminosity class lb 
jMafz-Apellaniz et alJ l2004f) , and as such is significantly 
cooler than the 04 prototype ^ Pup, which shows X-ray 
emission line profiles consistent with the wind-shock 
scenario. The wind mass-loss rate determinations for 
^ Ori are roughly a factor of two lower than those for 
^ Pup. Other important stellar and wind parameters 
taken from the literature are listed in Tabled The overall 
X-ray properties of ( Ori are quite typical of O stars 
(kTx < 1 keV, L^/Lboi ^ 10^1 iCassineUi fc Swank! 

| l983l: iBerghoefer. Schmitt. fc Cassinellil Il996l: 

IWaldron fc Cassinellill200ir 



3 PHENOMENOLOGICAL AND 

NON-PARAMETRIC ANALYSIS OF THE 
LINE PROFILES 

The simplest method, and a common mode, of examining 
hot-star emission line properties is the fitting of Gaussian 
line-profile models. These are convolved with the instru- 
mental response function and multiplied by the instrument 
effective area and fit to individual lines allowing for an as- 
sessment of the centroid shifts, line widths, and amplitudes. 
Indeed, this approach was taken in the paper that presented 
and fir st discussed the high-resolut io n Chandra spectra of 
C Ori JWaldron fc Cassinellil 1200 j) . IWaldron fc Cassinelli 
i2001I) reported significant broadening (velocity dispersion 
of 900 ± 200 km s~^), but also noted the generally symmet- 
ric appearance of the lines and reported a lack of any trend 
in line centroid shifts. 

We recapitulate this approach here, but also quantita- 
tively examine the quality of the Gaussian fits, including the 
distribution of the residuals. In Figure |21 we look at two of 
the strongest unblended lines in the spectrum, O vill Lya at 
18.969 A and Fexvil at 15.014 A. In these fits, shown in the 
top panels of each column, the centroid of the Gaussian was 
first fixed at the laboratory rest wavelength (the oscillator- 
strength-weighted mean of the two components of the Lya 
doublet in the case of the oxygen feature) and a power law 
was fit simultaneously to the weak contirmum. These fits are 
formally bad when analyzed using Monte Carlo simulations 
of the C statistic distribution - rejected at more than the 
90% level. There are clear indications of line profile asym- 
metries in the residuals of the Gaussian fit, in the sense one 
would expect from a wind-shock model, with a blue-shifted 
peak and steeper blue wings and shallower red wings. 

We next fit a Gaussian model with the centroid allowed 
to be a free parameter. This model (shown in the middle 
panel of each column in Figure|2l ^he line profile better, 
but there are clearly still systematic trends in the distribu- 
tion of fit residuals. Again, the actual line profiles have blue 
wings that are steeper than the Gaussians and red wings 
that are shallower. The Monte Carlo analysis of the C statis- 
tic distributions shows that these fits are better than those 
with the fixed Gaussian centroids, having rejection proba- 
bilities of only 68% and 73%, for the O vill and Fexvil lines, 
respectively. 

The widths and centroid shifts can be estimated from 
these Gaussian fits, even if the model is not ideal. For the 
oxygen Lya line, we find a best fit Gaussian (half width at 
half maximum) HWHM of 810 ± 30 km s"\ and a centroid 
blue shift of —150 ± 30 km s~^. Most other lines have even 
larger shifts, as can be seen in Table H in which we show the 
results of fits to seven emission lines in the spectrum. These 
values seem plausible in the context of the wind-shock sce- 
nario, although one might ask what values of the peak blue 
shifts and HWHMs would be expected. The estimated ter- 
minal velocity of the wind is, after all, twice the value of 
the derived HWHMs. The answer will depend on the spa- 
tial distribution of the X-ray emitting plasma, the veloc- 
ity distribution, and the degree of attenuation (see fig. 2 in 
lOwocki fc CohenI i200 J) !. We will show in the remainder of 
this section and the next one that there are quantitative 
indications of line asymmetries, even apart from the appli- 
cation of any specific wind model, and that an empirical 
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wind model does in fact fit tfie line profiles better than the 
shifted Gaussian model (the wind-profile model is shown in 
the bottom two panels of Figure |5] but not discussed until 
§4). Although we cannot reject the shifted Gaussian model 
with a high degree of certainty for any one emission line in 
the spectrum of C, Ori, the Gaussian fitting suggests some 
degree of line profile asymmetry and, more generally, that 
a more appropriate and physically meaningful model might 
improve the quality of the fits. 

But before fitting wind-profile models, let us first char- 
acterize the line profile shapes using a model-independent, 
non-parametric analysis. We do this by computing the first 
three moments of the observed line profiles, describing re- 
spectively the centroid shift, width, and asymmetry of the 
line profiles, as computed from: 
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Here a; is a dimensionless wavelength variable scaled to the 
terminal velocity of the wind, with the laboratory rest wave- 
length of each line set to a; = 0, as a; = (y ^'^'^ 

f{xi) is the number of counts in the i*'' bin of N total bins 
at scaled wavelength Xi. Note that we have not normalized 
the third moment in our definition, in order to make the 
calculation of its formal uncertainty more straightforward. 
The standard definition of the skewness, s, is related to our 
definition of the third moment according to 

Ma 



a^ie:=i/(^o 

We propagate the formal uncertainties for each calcu- 
lated moment from the Poisson errors on the total number of 
counts in each (scaled) wavelength bin. We note that we have 
not corrected for the instrumental broadening, which is quite 
symmetric and not very large compared to the observed 
line widths, and so will not affect the first and third mo- 
ments significantly. We also have not corrected for the weak 
continuum present under each line or for the wavelength- 
dependence of the detector effective area. But both of these 
factors are explored in quantitative detail in the next section, 
and are shown to be negligible. We list the values of the first 
and third moments for the stronger, unblended lines along 
with their formal uncertainties in Table |3 The second mo- 
ments are not listed, although they are quite large, because 
we have already determined from the Gaussian fitting that 
the lines are broad and in the moment analysis we cannot 
separate out the effects of physical broadening from instru- 
mental broadening. We use only the unblended lines in this 
analysis because the moment values have meaning only if 
they are calculated on a symmetric domain about a; = 0. In 
all cases we use the domain [—1 : 1] and assume a value of 
Voo = 1860 km s~^ for the wind terminal velocity. In Figure 
|S]we show two emission lines, with the moment- analysis do- 
mains indicated, along with the laboratory rest wavelengths 
and the values of the first moments. 



The numerical values of the first moments are straight- 
forward to interpret. They represent the position of each 
line centroid in units of x. The values of the third moments, 
however, are difficult to interpret by themselves. But their 
significance level in terms of formal uncertainties (i.e. their 
"sigma" levels, listed in the final column of Table I^J are 
the relevant quantity for assessing whether each line has a 
non-zero skewness (asymmetry) that is statistically signifi- 
cant. Unshifted and symmetric lines should have first and 
third moments that are consistent with zero. The emission 
lines analyzed for ( Ori are significantly blue shifted (neg- 
ative first moments), which is consistent with the results 
of the Gaussian fitting, but which contradicts the asser- 
tion of IWaldron fc Cassinellil (jio^) that there are no sys- 
tematic red shifts or blue shifts in the emission lines. It is 
also clear from the moment analysis that the lines are sig- 
nificantly redward skewed (positive third moments), gener- 
ally between the 1 and 2 sigma levels for each line. This 
asymmetry was not noted in the earl ier analysis, which re- 
lied o n "eyeballing" the Gaussian fits dWaldron fc Cassinelli 
120011) . A redward skewness (along with the blue shifted cen- 
troids) is exactly what is expected from continuum absorp- 
tion i n the context of a fast , spherically symmetric stellar 
wind JOwocki fc Cohenl200j) . The redward skewness comes 
about from the steep blue wing and the more extended, shal- 
lower red wing. 



4 WIND PROFILE MODEL FITS TO THE 
EMISSION LINES 

In the previous section we showed that there is evidence 
for blue shifting, redward skewness, and broadening in the 
X-ray emission lines of the O supergiant ( Ori. These re- 
sults are consistent with the expectations of a generic wind- 
shock picture. To augment this model-independent charac- 
terization of the net profile shift and skewness, and to de- 
rive physical information about the applicability of a wind- 
shock model, let us next fit a simple, empirical wind-shock 
line-profile model to the relatively strong lines in the MEG 
sp ectrum of C Qri. We use the empirical wind-profile model 
of lOwocki fc CohenI (1200 J) , which is physical, in the sense 
that it accounts for the Doppler shifted emission and radi- 
ation transport, including continuum attenuation, through 
a three-dimensional, spherically symmetric expanding wind. 
The parameters of the model have specific, physical mean- 
ings related to the spatial distribution of the hot plasma 
and the amount of absorption by the bulk, unshocked wind. 
The model is empirical, in that it does not posit any specific 
heating mechanism, and thus is applicable to a wide range of 
possible wind-shock (and even coronal) scenarios for X-ray 
emission. 

The goal of fitting the wind-profile model is thus to 
constrain the physical parameters of the wind emission and 
absorption for each strong emission line in the Chandra spec- 
trum of Ori. Theorists may then compare the predictions 
of any number of specific models or numerical simulations 
to the physical parameter values we derive. Furthermore, 
our fitting of wind-profile models allows us to quantify the 
amount of asymmetry in the line profiles and relate the 
asymmetry, quantitatively, to the amount of wind absorp- 
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tion, through the optical depth parameter of the wind- profile 
model, 

nil 

where k is the absorption opacity, and M is the mass-loss 
rate. Note that bound-free continuum absorption is the dom- 
inant opacity source. Physically, r* represents the optical 
depth along a central ray from infinity to the stellar sur- 
face radius R,, in the simplified case that the wind veloc- 
ity is constant at the terminal value, v^a- In this simplified, 
constant-velocity value of r* > 1 also represents the 

radius of unit optical depth, R\, expressed in units of R*. 
The wind-profile model assumes that the hot plasma is dis- 
tributed throughout the wind, above some minimum radius, 
-Rmin, and that its filling factor is proportional to the am- 
bient wind density multiplied by an additional power-law 
factor, / (X r^' (thus falling off as l/iir'^^''). The other 
two interesting parameters of the model are thus _Rniin/R* 
(sometimes expressed as Umax = R./Rmin) and q. The nor- 
malization of the profile is the fourth, and final, parameter. 
There is an implicit assumption that there are enough differ- 
ent regions of hot plasma that the wind can be treated as a 
two-component fluid, comprising a bulk, cool (T « Tch), X- 
ray absorbing component, with a hot. X-ray emitting com- 
ponent smoothly mixed in. The minimum radius of the hot 
plasma distribution is motivated by numerical simulations 
that show that large shocks tend not to form until the wind 
flow has reached at least several tenths of a stellar rad ius 
JCohen et al.lll99el : ICooperlll995 ^^^^^^^j^W^L 

A s discussed in further detail ii nOwocki^^Cohenl 
l|200ll) . the line profile is computed from the integral 

y.-(<7+2) 

^"""y (I -R^/r)^/^ exp[-r(M.,r)]dr, 

where = max[iimin, R./(l - Ixl^/")], fi^ =x/{l- R*/v f, 
and T{^,r) (which is proportional to r*) is the optical 
depth along the observer's line of sight at direction co- 
sine fj, and radial coordinate r. Here the scaled wavelength, 
X = (A/Ao — l){c/voa), is the same quantity we used in the 
moment analysis. The parameter (3 is the usual wind accel- 
eration parameter, from v = Voo(l — R./r)'^. The governing 
equation for Lx must be solved numerically for all /3 ^ 0. We 
set /3 = 1 in all of our fits. We include a power law contin- 
uum model in all the fits we performed in XSPEC. Finally, 
we note that this profile model implicitly assumes spherical 
symmetry and a smooth wind flow. 

Again, this wind-proflle model is both physically mean- 
ingful and widely applicable to a range of different physical 
models of X-ray production, inc luding coronal models (see 
fig. 2 in lOwocki fc CohenI (|2003) for a graphical exploration 
of the effects of choosing different model parameter values 
on the line profile shapes, and fig. 4 in the same paper for a 
comparison of wind-shock and coronal model parameters). 
The larger -Rmin is and the smaller q is, the broader the line 
profiles tend to be. We note that for a wide range of realistic 
choices of these parameters, the characteristic width of the 
resulting profiles is equivalent to roughly half the terminal 
velocity, consistent with the half- widths we derived from the 
Gaussian fits in the previous section. Increasing the wind op- 
tical depth parameter, r«, tends to make the profiles more 



narrow, more blue shifted, and more asymmetric. A model 
with a relatively small Rmin value and a negligible r« pro- 
duces a profile that is similar in shape to a Gaussian. 

We fit this wind-profile model to each strong line in 
the Ori MEG spectrum, allowing all four adjustable pa- 
rameters (r«, -Rmin, q, and the normalization) to be free, 
in conjunction with a power-law component to model the 
weak continuum emission. For several line complexes, we 
fit multiple profiles simultaneously to account for blending. 
This included the helium-like resonance and intercombina- 
tion lines of oxygen. We do not give flts to the other helium- 
like complexes in the ^ Ori Chandra spectrum, as flts to 
these complexes are reported elsewhere iLeutcriea'ge r et alJ 
l200d) . For the Fexvil lines at 17.051 A and 17.096 A (3G 
and M2, respectively), which we flt simultaneously, because 
they are quite blended, we fixed the relative normalizations 
to I(M2)/I(3 G) = 0.8, consistent with HULLAC calcula- 
tions (Maucho . Liedahl. fc Fournicr 2001) and with the val- 
ues generally observed in stars. When fitting these lines, and 
also other line blends, we tied the three primary parameters 
of the line-proflle model - r«, -Rmin, q - together. Ultimately, 
we report here on the flts to the nine lines in the spectrum 
that provide meaningful constraints to the model parame- 
ters. Note that these are not the same set of lines to which 
we fit Gaussian models, as it was easier to get meaningful 
Gaussian flts to several weaker lines and it was harder to 
get meaningful Gaussian fits to the blended lines. 

As mentioned previously, we first carried out this mod- 
eling using the same procedure, implemented in Mathe- 
matica, that we em ployed in our earlier analysis of ( Pup 
jKramer et alJbOOaTl . We then repeated the modeling using 
a custom-written module in XSPEC, which allowed us to 
include a continuum emission component in the modeling 
and use the exact instrumental responses. This also enabled 
us to simultaneously fit multiple models to line blends. The 
two methods gave very similar results, and the XSPEC flts 
were insensitive to both the choice of continuum model and 
the wavelength range included in the flt. The results of the 
XSPEC model-fltting are summarized in Table 2] and the 
best-flt models, superimposed on the data, are shown in Fig- 
ure 2] and the bottom row of Figure |5| for the nine lines, in 
seven complexes, for which we could obtain meaningful fits. 
The wind-profile model does indeed provide better flts to the 
stronger, unblended lines than does the Gaussian model, ac- 
cording to the Monte Carlo simulations of the distribution 
of C statistic values. The goodness of fit values (expressed as 
a percentage of the Monte Carlo simulations that gave a C 
statistic as good as or better than that derived from the fit 
to the actual data; lower percentages are better) are listed 
in Table |1| All the wind-proflle flts are formally good. 

We calculated errors on the derived model parameters 
by using a three-dimensional grid of models in the parameter 
space of interest (r*- -Rmin - q) and applying a AC criterion 
appropriate for jointly-distributed uncertainties for three pa- 
rameters, and report the maximum extent of this confidence 
region in each of the three parameters as the formal un- 
certainties on the derived parameters. These are the values 
listed in Table |3] (90% confidence limits for one parameter 
of interest - AC = 2.71), and shown, for two particular fits, 
in Figure]^ 
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5 DISCUSSION 

We summarize the derived model parameters and their un- 
certainties for each line in Figure El This figure shows that 
there are no strong trends in any of the wind-profile model 
parameters with wavelength (or any other characteristic) of 
the emission lines. Fitting a function linear in wavelength 
to the uncertainty-weighted model parameters shows con- 
sistency with a constant function (at the 95% confidence 
level) for each of the three parameters. Note that the r, 
point for the O vil complex near 22 A must be excluded for 
this statement to be true. We discuss this outlier in terms 
of the wavelength dependence of the wind opacity near the 
end of this section. 

The fitting results shown in Figure |S| present a consis- 
tent picture of a line profile model with t, ~ 0.25 to 0.5, 
an onset radius, i?min ~ 1.5R, , and a constant filling factor 
(g « 0). These are all reasonable parameters in the context 
of the general instability-driven wind-shock model, though 
the r* values are small compared to the expectations of wind 
theory, which we elaborate on below. Finally, we note that 
most of the lines cannot be well fit by models with no wind 
absorption (r, = is ruled out), which is consistent with 
the inability of Gaussian models to provide good fits and 
also with the non-zero third moments of the line profiles, as 
discussed earlier. Looking at the situation from a different 
point of view, upper limits on the wind absorption are above 
r* ~ 0.5 for all but one line complex in the spectrum. The 
unmistakable conclusion is that the Chandra spectrum of 
C, Ori is consistent with a moderate amount of wind absorp- 
tion (as well as the expected degree of broadening from an 
embedded wind source), and that at least some wind atten- 
uation is demanded by the data. 

The derived -Rmin and q values are consistent with 
the numerical simulations of the line-force instability wind 
shocks, inferred from simula tion output shown in various fig- 
ures inlCoopeJ lll99 6l: Co hen et al.Nl99ell : iFeldmeier et all 
^1997b^ : IOwocki fc RunacresI ll200al . These trends are also 
qualitatively understood from a theoretical point of view. 
The strong, relatively symmetric diffuse (scattered) radi- 
ation field near the photosphere inhibits the line-force in- 
stability and thus the formation of strong shocks near the 
photosphere, and the filling factor is not strongly depen- 
dent on radius because although the propensity of shocks to 
form eventually falls off with distance from the photosphere, 
the cooling timescale for shock-heated plasma increases with 
distance. 

Given the spatial distribution of hot plasma derived 
from the line-profile fits, the continuum attenuation by the 
overlying cool wind is governed by the mass-loss rate and 
wind opacity. In the model we have employed, the overall 
wind attenuation is characterized by the optical depth pa- 
rameter, 

47ruooRt 

Using mean values from Table Q and a wind opacity value 
of K « 125 cm'^ g~^, we expect r« ~ 3. The va l ue for the 
wind opacity is taken from fig. 4 in lCohen et alj il99(tl . Of 
all the values that go into this calculation, the mass-loss rate 
is probably the most uncertain, followed by the wind opac- 
ity and the star's radius. The terminal velocity is probably 



known to within ten or twenty percent (which is the range 
of values found in the literature). 

Thus, the value of the wind optical depth parametrized 
by r», as derived from the observed X-ray line profiles, is 
about an order of magnitude lower tha n the expected value . 
This is similar to what is seen in C Pup jKramer et all2003ri . 
where the observed value of r, is almost an order of mag- 
nitude lower than expected (there is also a fair amount of 
uncertainty in the relevant properties of C, Pup). The ex- 
pected r* value for C, Pup is about a factor of two larger 
than that for C, Ori, primarily because of the earlier type 
star's larger mass-loss rate. 

The fact that the X-ray line profiles of C, Pup, and now 
C, Ori, indicate lower than expected wind optical depths 
is consistent with recent work that suggests that O star 
mass-loss rates may have been overestimated by a factor 
of three or more, an d perhaps up to an order of magnitude 
iBouret et al.lboo'sh due to clumping (which affects density 
squared mass- loss diagnostics, such as radio free- free and 
Ha emission). This result is not inconsistent with the tra- 
ditional UV absorption-line based mass-loss rate estimates 
of hot-star winds, which have always been subject to uncer- 
tainty due to the difficulty of reliably accounting for ioniza- 
tion distribution effects. In fact, other recent work, focusing 
on far-UV absorption line studies as diagnostics of mass loss 
and wind ionization in many O and B supergiants, indicates 
that mass-loss rates based on UV absorption line analysis 
may be overestimated by as much as an order of magnitude 
jFullerton. Massa. fc Priniall2006ll . 

Recent detailed multi-wavelength modeling of a large 
sample of O giants and supergiants (but not including Ori) 
indicates mass-loss rate overestimates of at least a factor of 
two, assuming that the far wind, where the radio free-free 
emission arises, is unclumped, and more than a factor of two 
if the far wind is significantly clumped 4Puls et al. 200t|). 
This work also shows that there are star-to-star variations 
in clumping factors and a somewhat strong radial depen- 
dence - at least for some stars - to the clumping factor. We 
note in this context that the X-ray profiles will in general 
be most sensitive to clumping in the region near and just 
a bove -Rmin (so, R ~ 2Rt , which is "region 2" in the analysis 
of iPuls et ail (120061) ). These various threads of evidence for 
lower mass-loss rates are also consistent with the energy bud- 
get analysis of wind-b l own bubbles a nd superbubbles (see, 
e.g.. lNaze et al.l feOOSH : ICooper et alJ 12004.1 and references 
therein; but see also iFrever^Henal^r^^^A^ ll2006ll for the 
role played by the swept-up wind from earlier evolutionary 
stages in Wolf-Rayet bubbles). 

Even apart from its effect on mass- loss rate estimates, 
clumping itself has the potential to reduce the effective 
opaci t y of a stellar win djjFel^ieigr^Osldnflva^ fc HamannI 
120031: lOskinova. Feldmeier. fc HamannI |2004^ . This ef- 
fect might more accurately be termed "porosity," as 
it presumes the existence of a low-density interclump 
channels that can pote ntially allow photons to escap e 
the wind more easily jOwocki, G^de^fc Sh avivl l2004l) . 
lOskinova. Feldmeier. fc Ifamam!n2005n2006t) have recentlv 
computed X-ray line profiles for specific models incorporat- 
ing geometrically thin, radially compressed shells, compar- 
ing their results with observed profiles for several hot stars 
(including C, Ori). Within the assumptions in their model, 
these authors show that a very optically thick emission line. 
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with a significant skewness in a smooth wind, can be made 
moderately more symmetric with an interclump spacing of 
0.2 R,, and can be made nearly s ymmetric w ith an inter- 
clump spacing of 2 R« (see fig. 1 in iiOskinova et al. (2005)). 
Similarly, using a parametrized model of isotropic clump- 
in g, based generally on the porosity formalism introduced 
bv lOwocki et al.l J2004) . lOrocki fc C ohen (200^ find that 
obtaining symmetric X-ray emission profiles from an oth- 
erwise optically thick wind requires a quite large "porosity 
length" h = £/ f, where £ represents the characteristic clump 
scale, and / if the clump volume filling factor. Specifically 
it requires h of order the local radius r. 

In this context, it is thus important to stress that while 
the mass-loss overestimates due to clumping depend only on 
the density contrast between the clumps and the interclump 
medium (and thus the volume filling factor) , for porosity to 
affect the line profiles directly, the density contrast must be 
accompanied by a large clump scale, or interclump spacing. 
The most sophisticated numerical treatment of the line-force 
instability shows structure on small {£ ^ R») spatial scales, 
with only moderately compressed volume filling factors (/ ~ 
0.1 iPcssart fc Qwocki .2Q0a) ). 

These results make it difficult to see how the wind inho- 
mogeneities produced by the instability and which, presum- 
ably, are directly related to the shock-heating responsible for 
the X-ray emission itself, could lead to a significant porosity 
effect on the X-ray line profiles. In light of the several in- 
dependent lines of evidence for lower O star wind mass-loss 
rates, we suspect that lower wind column densities are the 
cause of the order of magnitude discrepancy between the r* 
values we derive in this paper from fits to the emission lines 
in the C handra spectrum of C Ori and the similar results de- 
rived bv lKramer et al.) (|2003) for ^ Pup. Clumping and the 
associated porosity may play some role, but for that role 
to be significant, the porosity length in O star winds must 
be large - of order the local radius - and the combination 
of these two effects must reduce the effective wind optical 
depth by an order of magnitude. 

The results from our X-ray emission line profile analy- 
sis should be consistent with other aspects of the Chandra 
observations. The emission measure an d temperature infor- 
mation derived from the ob servatio ns (|CassineULfc_Swanlj 
Il983l: iBerghoefer et alJ IT996: Wald ron fc Cassinellil l200l 
are typical for O supergiants and do not provide 
any significant constraints on the interpretation of the 
line profiles, aside from simply being broadly consis- 
tent with the expectations of the standard wind-shock 
scenario. The most constraining specific X-ray diagnos- 
tic in conjunction with the emission line profiles is the 
forbidden-to-intercombination emission line ratio in the 



helium-like isoelectronic sequence 


Gabriel & Jordanlll969l 


iBlumenthal. Drake. & Tuckeil Il972 


). In the presence of a 



strong UV field which can drive photoexcitation of electrons 
from the upper level of the forbidden line to the upper level 
of the intercombination line (2s ^Si - 2p "^^1,2) and thus 
reduce the f/i line ratio, it can be used as a diagnostic of 
the UV mean intensity and thus of the distance of the X-ray 
emitting plasma from the photosphere. 

The initial work on the several helium-like f/i ratios 
seen in the Chandra spectra from 1^ Ori showed that most 
of the helium-like ions were far from the photosphere, con- 
sistent with those ions being embedded in the stellar wind. 



but that the Sixill f/i ratio implied a location only slightly 
above the photosphere, which would generally be considered 
too close to t he star to b e consistent with any wind-shock 
scenario ( Wa ldron fc Cassinelli 2001') . However, a recent re- 
analysis of these same data showed that all the helium like 
ions, including Sixill, a re consistent with an on set radius 
{Rmin) of about 1.5R, iLeutenegger et al]l2006l) . This re- 
sult is, of course, completely consistent with those we report 
here for the emission line profiles of nine other lines in the 
Chandra data. 

We can also consider trends in the derived wind profile 
parameters within our dataset. One might expect different 
lines to have different morphologies and thus different model 
parameters either because different ions form at different 
temperatures and thus sample different shocked regions or 
because lines at different wavelengths have differing amounts 
of wind attenuation due to the wavelength dependence of 
the opacity of the bulk, cold wind. Regarding the first pos- 
sible effect, we note that numerical simulations show a rel- 
atively constant rms velocity dispersion with ra dius, once 
sh ocks begin to form ('Runa cres fc Owockill2o61h . Figure 5 
in lRunacres fc Owocki. (.20Q2il shows, in detail, a very rapid 
rise in the velocity dispersion, followed by a very shallow 
fall-off with radius. 

Regarding the second effect; that of wavelength- 
dependent attenuation, photoionization cross sections of cos- 
mically abundant plasma do have a strong wavelength de- 
pendence over a large range of wavelengths. However, this 
effect is more complex when the plasma is ionized, as it is 
even in the "cold" component of a hot-star wind. Further- 
more, the lines we analyze in this paper span only a factor 
of t wo in wavelength. L ooking at the wind opacity in fig. 
4 of ICohen et"ai] (Il996h . we can see that the values of the 
wind opacity range only over about a factor of 2 from 600 
eV (roughly the photon energy of the O vil lines near 22 A 
which are the longest-wavelength lines to which we fit the 
wind profile model) to 1000 eV (roughly the photon energy 
of the Nex Lya line, which, at A = 12.134 A, is the short- 
est wavelength line we discuss here). The variations in the 
wind opacity on this relatively small wavelength range are 
complex and not monotonic because of the dominance of 
photoionization edges of oxygen (O"*"^ through O"""^). The 
appearance of these edges breaks up the usual fall off 
in opacity, and over this relatively small wavelength range, 
makes the opacity roughly constant. If anything, the longest 
wavelength lines in our data (the O vil lines near 22 A) are 
subject to less attenuation than the shorter wavelength lines, 
by virtue of their being longward of the oxygen K-shell edges 
(and, in fact, this emission feature has the lowest upper limit 
to the r, parameter of any of the lines we fit). In any case, 
there are no significant systematic trends in any of the three 
wind profile model parameters. As we discussed above, a sin- 
gle value of each parameter is consistent with all the data. 
So, although higher signal-to-noise data in the future may 
reveal a significant trend, no ne is seen in these data. We 
should point out, though, that lOskinova et~ al. (2006) noted 
that radiation transport through a medium with completely 
optically thick clumps will not only reduce the effective wind 
opacity, but will make the opacity effectively gray. Interpret- 
ing the wavelength dependence of line profile morphologies 
- or lack thereof - however, requires both a detailed eval- 
uation of the wavelength-dependent atomic opacity and its 
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uncertainty, and also statistical fitting of whatever line pro- 
file model may be appropriate along with formal constraints 
on confidence limits of the parameters of that model. 

Finally, we note that each line or line complex is well 
fit, in a statistical sense, by the relatively simple, spheri- 
cally symmetric wind-profile model we employ here. Future 
higher resolution and/or higher signal-to-noise spectra could 
show evidence for signatures of wind asymmetry or of time 
variability in the line profiles (perhaps much like DACs see 
in UV absorption lines from the winds of hot stars or like 
moving emission bumps seen in WR spectra). There is, how- 
ever, no need at this point to invoke either of these effects 
nor any others that go beyond the basic model we have used 
here. 



6 CONCLUSIONS 

The fundamental observational conclusions of this work are 
that the X-ray emission lines of the late O supergiant C, Ori 
are broad, blue shifted, and modestly asymmetric, which is 
qualitatively consistent with the general picture of hot, X- 
ray emitting plasma embedded in an expanding, spherically 
symmetric stellar wind. These results come both from fit- 
ting a physics-based empirical wind-profile model to nine 
emission lines in the Chandra MEG spectrum, and also 
from attempts to fit Gaussian line-profile models and a non- 
parametric analysis of the line shapes via the calculation of 
the first three moments of seven unblended lines. 

There is no need, based on the observed line profiles, to 
invoke ad hoc coronal emission or other non-standard X-ray 
production mechanisms. However, the amount of attenua- 
tion by the bulk, cold stellar wind is significantly less than 
would be expected by a simple application of the assumed 
mass-loss rate, standard warm plasma opacities, and the as- 
sumption of a spherically symmetric, smooth stellar wind. 
Qualitatively, this result is consistent with the results of a 
similar analysis of the Chan dra spectrum of the early O star, 
CPup (E ramer et alJbOOaTl . And the smaller-than-expected 
wind attenuation leaves an observational signature that ex- 
plains why previous studies, in which Gaussian profiles were 
fit and then analyzed "by eye," did not identify the signature 
of wind attenuation. The emission lines, though significantly 
blue shifted, are only modestly asymmetric, and in fact, any 
individual line can be at least marginally fit by a blue shifted 
Gaussian. For the strongest lines, however, there is a sig- 
nificant improvement in the fits based on the wind-profile 
models as compared to those based on Gaussians. 

These results, taken together with the earlier ones on 
the X-ray line profiles of Pup, indicate then that the 
standard wind-shock scenario is adequate for explaining the 
high-resolution X-ray spectra for normal O supergiants. Un- 
usual hot stars, such as 9^ Ori C and r Sco, do not fit 
into this paradigm, perhaps because of their extreme youth, 
but there is no reason, especially now t hat the helium-like 
f li li ne ratios have also been reanalyzed iLeutenegger et all 
l2006f) . to suppose that all hot stars, with the sole exception 
of C Pup, pose an insurmountable challenge to the wind- 
shock model of X-ray production. That being said, the wind- 
shock model still has various difficulties in accounting in de- 
tail for the observed trends in X-ray properties among OB 
stars, and there are many open questions about the specific 



ingredients of a correct wind-shock model. But the nature 
of X-ray emission line profiles in O supergiants, while pro- 
viding some interesting constraints and presenting a puzzle 
about wind optical depths, does not require us to completely 
discard the wind-shock paradigm or lead us to invoke coro- 
nal models for explaining hot-star X-ray emission. The lower 
than expected wind optical depths derived from the X-ray 
line profiles do, however, add to the debate about O star 
mass-loss rates and the role of wind clumping. 
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Table 1. Stellar properties of C Orionis from the literature. 
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2100 


Voels et al. (1989) 


34 


24 


— -0.27 






Wilson and Dopita (1985) 


25 


20 




1.58 


2190 



Table 2. Gaussian line profile fits to the emission lines. 



Ion Ao Ccntroid HWHM 

(A) (km s-i) (km s-i) 



Nvii 


24.781 


-110 ± 140 


1380111° 


O VIII 


18.969 


-150 ±30 


810 ± 30 


Ovii 


18.627 




500lif 


O viii 


16.006 


-ioo«o 


880 ± 70 


Fe xvii 


ir5.oi4 


830l|0 


Nex 


12.134 


-150 ±50 


980 ± 50 


Ne IX 


11.544 


-390 ± 140 


1360ll«° 



Table 3. First and third moments of the emission line profiles. 



Ion 


Ao (A) 


Ml 


Ml/uncert. 


M3 


M3/uncert. 


O VIII 


18.969 


-0.0818 ± 0.0135 


-6.08 


7.9914 ± 4.8511 


1.65 


Ovii 


18.627 


-0.1138 ±0.0416 


-2.74 


4.5943 ± 2.9983 


1.53 


Fe XVII 


16.780 


-0.1652 ±0.0229 


-7.32 


7.5967 ± 3.9259 


1.94 


O VIII 


16.006 


-0.0464 ± 0.0247 


-1.88 


3.9309 ±3.8114 


1.03 


Fc XVII 


15.014 


-0.0792 ± 0.0173 


-4.58 


12.5198 ± 5.5363 


2.26 


Nex 


12.134 


-0.0801 ± .0194 


-4.13 


10.1529 ± 5.4936 


1.85 


Ncix 


11.544 


-0.1108 ±0.0368 


-3.01 


2.8778 ±4.0377 


0.71 



Note: Ml and M3 arc the first and third moments of the line profiles, respectively. The following columns show the ratio 
of the values of these moments, for the indicated unblended lines, to their formal uncertainties. We interpret the values in 
these columns as significance indicators of the first and third moments' deviation from zero, as described in Sec. 3. 



Table 4. Wind profile model parameters fit to the data. 



Ion 


Ao (A) 


1 


R-miii/ R* 




Goodness of fit' 


Ovii 
O VIII 
Ovii 
Fe XVII 
O VIII 
Fe XVII 
Nex 


21.804, 21.602 

18.969 

18.627 
17.051, 17.096 

16.006 

15.014 

12.134 


-0.30l;2^ 
-0.12l:g 
0.39+^7^ 

-o.4ii:f, 
-o.4ii:« 

-0.47l-jg 
-0.50i;»5 


I.66+I3 

1 61+-" 

l.ol_ j2 

1.291:29 

1.28+13 

1-511-25 
1.37115 

1.55i;^? 


o.o6l:J^ 

0.26+-2'] 
1 34+1.86 

0.76+H 
0.27+-« 
0.58+-^; 
0.45i;« 


0.33 
0.67 
0.49 
0.40 
0.53 
0.29 
0.14 



"Fraction of Monte Carlo simulated datasets that gave a C statistic as good or better than that given by the best-fit model 
and the data. This can be interpreted as a rejection probability. Lower values indicate better fits. 
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Figure 1. The 



MEG spectrum of C, Ori, witli negative and positive first orders from both observations (Obs. IDs 610 and 1524) 



coadded. 
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Figure 2. Best-fit models superimposed on the observed Ovill Lya line (left-hand column) and the Fexvil 15.014 A line (right-hand 
column). Error bars here and in other figures are calculated from the total source counts per bin, assuming Poisson errors. The fits shown 
in the top row are for a Gaussian model with the line center fixed at the laboratory rest wavelength. This fits shown in the middle row 
are for the Gaussian model with the centroid treated as a free parameter. The Gaussian fits in the first two rows are discussed in Sec. 3. 
The fits shown in the lower panel are for the wind-profile model discussed in Sec. 4. 
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Figure 3. The same two representative emission lines shown in Figure|3] with their centroids as determined from the moment analysis. 
Ovm Lyman-alpha at 18.969 A (left) has a centroid (first moment) 6 sigma from the laboratory rest wavelength, and a positive third 
moment (red skewed) that is significant at the 1.7 sigma level. The Fexvil line at 15.014 A (right) has a significantly negative first 
moment (5 sigma) and a third moment that is positive at the 2.3 sigma level (see Table IS! . In both panels, the solid vertical line is the 
laboratory rest wavelength, while the dashed line to its immediate left represents the first moment. The other two dashed lines represent 
the blue and red limits over which the moment analysis was performed (a; = — 1, 1). 
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Figure 4. Best-fit wind-profile models for five lines (or line complexes): Nex 12.134 A, Ovili 16.006 A, Fexvii 17.051 A and 17.096 A, 
Ovii 18.627 A, and Ovil 21.602 A and 21.804 A. 
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Figure 5. The 68% and 90% confidence regions in the parameter space of the wind-profile model, for the two lines shown in Figure|21 the 
Ovill Lya line (left-hand column) and the Foxvil 15.014 A line (right-hand column). The best-fit model parameters are indicated by the 
asterisks. Note the correlation between q and iJmin ("max)- For each 2-D slice of parameter space shown here, the other model parameters 
are optimized (i.e. free) while models are fit for a grid of the two displayed parameter values. The contour levels thus correspond to AC 
values appropriate for two parameters of interest: 2.30 and 4.61. 
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Figure 6. The best-fit model parameters - r, , q, and -Rmin - for each line complex we fit with a wind-profile model. The error bars 
represent the 90% confidence limits on the one parameter of interest in each panel (AC = 2.71). 



